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We propose an optomechanical setup where the activation of entanglement through the pre-availability 
of non-classical correlations can be demonstrated experimentally. We analyse the conditions under which 
£SJ . the scheme is successful and relate them to the current experimental state of the art. The successful 

activation of entanglement embodies an interesting alternative to current settings for the revelation of 
fully mechanical nonclassicality. 
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r Nderstanding quantum correlations and their role is one of the key goals of the current research in modern quantum 
' physics (l]] . Although much of the attention has so far been focused on entaglement as the form of quantum correlations to 
look at when dealing, for instance, with quantum computational speed-up and better-than-classical communication capabilities, 
such a perspective has been challenged by the identification of other forms of non-classical correlations going beyond entangle- 
ment |2|47l]. Needless to say, as entanglement is the unique form of quantum correlations for pure states, the arena were the role 
(if any) that such quantum correlation should be investigated is the one embodied by mixed quantum states. For such states, 
entangled is no longer synonymous of non-classical and other quantifiers have been proposed, each striving to grasp specific 
Q-* aspect through which quantumness of correlations manifests itself. Among them, quantum discord [2] has so far enjoyed a 
l— ~ 1 ■ growing popularity, notwithstanding the difficulties inherent in its analytic formulation even for simple two-qubit states, due to 
the intriguing im plic ations it is alleged to have for the speed-up of some protocols for quantum computing [8] and its operational 
interpretation ffifioll. 

' Remarkably, it has been very recently realized that an effective way to understand discord and its role in quantum information 
f-*- , processing could be the establishment of tight relations with entanglement itself. In particular, Refs. Ill ll Il2ll have shown that, 
^sO 1 in the n-qudit scenario aided by unitary quantum gates, discord might be interpreted as the resource whose availability enables 
physically relevant effects, such as the production of an entangled state useful for quantum information processing, otherwise 
prevented when only classical correlations are at hand. Besides providing an alternative definition of classically correlated states, 
t-H | these results give a direct link between broader quantum correlations and the resource for quantum information processing. To 
date, although discord has been successfully extended to the domain of Gaussian continuous variable states lfl3l [Till , it is not 
\ I ' known whether a similar schemes hold in the continuous variable arena as well. 

, Here we provide strong evidence that, indeed, an entanglement activation protocol based on the pre-availability of discord 
can be formulated for registers and ancillae having continuous spectra. In order to do so, we adapt the original formulation of 
Ref. lfl2Tl to a realistic system of cavity quantum optomechanics that includes, ab initio, environmental losses and decoherence, 
thus extending the scheme itself from the unitary to the open-system domain. We use two independent optomechanical cavi- 
ties B15I417I1 : the register is embodied by the mechanical mirrors, the ancillae by the cavity fields. We show that pre-avaliable 
quantum correlations in the state of the mechanical subsystem can always be transformed into entanglement involving the 
mechanical and the optical systems. Besides its fundamental relevance, our study can be seen from the perspective of effective 
diagnostics on the purely mechanical state: the current experimental progresses are such that optomechanical entanglement 
will soon be within reach thus proving the potential for experimental quantum information science held by the optomechanical 
arena lfl8Tl2lll . By revealing the appearance of optomechanical entanglement, our scheme will give us indirect evidence of 
(much harder to access) nonclassicality at the mechanical level. Therefore, this study contributes to the grasping for a full 
comprehension of the interplay and relation between discord and entanglement and proves the possibility of activating the 
latter through the former in an experimentally realistic CV context. Furthermore, it opens up the way to further uses of 
optomechanical systems as nodes in quantum networks. 

Results 

Entanglement activation scheme. For clarity, here we briefly summarise the working principle of the entanglement activation 
scheme put forward in Ref. Il2tl and illustrated in Fig. 1 (a). The protocol is performed between two registers, the system 
register S (composed by Sj with j = 1, n) and the ancilla register A (with elements A,). The two registers are initialised 
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FIG. 1: Schemes of principle of the activation protocol, (a) The register composed of systems Sj (j = 1, 2) is initialized in a separable 
state and subjected to a set of local rotations (represented by the circled arrows). System Sj then evolves jointly with the ancilla Aj according 
to the completely positive map $. If the initial state of the register has non-zero discord, the output state after the application of the maps is 
entangled with respect to the system-ancilla bipartition (Si&X^li^b). The arrow of time is also shown, while D stands for the presence of 
discord, (b) Each system-ancilla unit is embodied by an optomechanical system consisting of an open optical resonator with a light, movable 
end mirror, pumped by an intense light field. The field-mirror radiation-pressure coupling, together with the lossy mechanisms affecting the 
device (see body of the manuscript for details) realize the map <E>. Separable yet discorded initial mechanical states are prepared by driving the 
cavities with a two-mode squeezed vacuum of proper squeezing parameter. 

in the factorized state ps ® pA, where ps is a separable state of the system register and pa— ®|Li Pa * s tne tensor product 
of n identical states of the ancilla register. A demon applies local random unitaries Ui on each Sj, which then undergo a joint 
evolution with the homonymous element of the ancillary register Aj. In general, such evolution is described by a completely 
positive dynamical map The key statement of reference 11211 is that iff ps is classically correlated (i.e., having zero discord) it 
is always possible for the demon to find a configuration Ud of demon's unitaries such that (S>^ =1 ^i[Ud(ps ® Pa)U\] is separable 
across the system-ancilla split. The result has been found under the assumptions of i) registers composed of qubits, ii) a rigid set 
of unitary gates applied to the Sj-Aj systems and Hi) ancillae initially prepared in a fiducial reference state. 

Activation scheme in a realistic optomechanical device. We turn the key claim in Ref. {T20 somewhat upside-down and 
provide the following 

Statement: Consider a register composed ofn non-interacting optomechanical systems, each interacting with a CV ancilla pre- 
pared in a fiducial initial state. If the elements of the register initially share non-classical correlations, dynamical entanglement 
in the register-ancillae split can be created robustly against adversary effects of a demon that acts locally on the state of each 
element of the register itself. 

As anticipated, the physical system comprises two optical resonators, each endowed with a mechanical end-mirror and inter- 
acting via radiation pressure with two optical modes IU5I - U7I1 . Each cavity is pumped by an external laser field that is quasi- 
resonant with the cavity frequency uj % c . The Hamiltonian of each optomechanical device, in a frame rotating at the frequency ujl 
of the external lasers (assumed to be the same), reads 



where (pi) is the position (momentum) quadrature of the i mechanical system, Cj (cj) is the annihilation (creation) operator 
of the i th cavity field (whose photon-number operator and energy decay rate are hi and m, respectively), and S 1 —uj 1 c ~ lol is the 
pump-cavity detuning. Each cavity has length Li, so that Xi = ^ /-^* * s me associated cavity-mirror radiation pressure coupling 
rate. The frequency and mass of the i th mechanical oscillator are uj l m and m*. Finally, by calling Vi,2 the pumping power of 
each driving laser, we have E.- L — yj2n.iPi/ (hij\f). 

The mechanical system is affected by decoherence induced by the thermal Brownian motion of the mechanical oscillators, 
which are damped at a rate 7^. Together with the cavity losses that have already been mentioned, this makes our system 
explicitly open and the dynamics undertaken by the joint optomechanical device properly described in terms of the completely 
positive map p out — <E>i = i^i[pin\ with pi n the initial state of the joint mirror-light system. We provide details on the explicit 
form of each <&i in the Methods section. 

Activating entanglement from quantum discord. We now show how quantumness of correlations between the mirrors guar- 
antees the generation of dynamical entanglement in the mirrors-light bipartition. In the Methods section we recall the definitions 
of discord, quantifying quantumness of correlation, and entanglement for a two-mode Gaussian CV states. Before discussing 
the actual results, we make a preliminary analysis of the entanglement created in a single optomechanical cavity between the 
mirror and the cavity mode as a function of the temperature of the Brownian bath the mirror is interacting with. 
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FIG. 2: High temperature entanglement activation, (a) The maximum value of dynamical entanglement created over time between the 
mechanical and optical mode as a function of the temperature of the brownian bath, (b) Optomechanical entanglement created between the 
field of one cavity and the corresponding mechanical mirror as a function of the interaction time. The entanglement in the other optomechanical 
subsystem shows an identical behavior. The initial state is the tensor product of coherent states of the cavity fields and a separable yet discorded 
state of the mirrors obtained using the protocol in Ref. 12911 . The inset shows the distribution of output entanglement achieved using 10 3 
randomly rotated initial covariance matrices vrjo corresponding to the same parameters used in the body of the manuscript. In all panels, 
A = LJ m . 

Radiation pressure is responsible for an in principle entangling dynamics between the mechanical mode and the cavity 
field I22H28I1 . However, such entanglement is the result of a delicate trade-off between the strength of the optomechanical 
coupling, the mechanical mass, the mechanical/optical quality factor and, quite crucially, temperature. The generation of en- 
tanglement during the transient evolution does not guarantee its persistence at the steady state achieved by taking the limit for 
t — > oo in Eq. In particular, while the steady-state entanglement will depend only on the relative ratios of the critical pa- 
rameters identified above, the dynamical one will have a strong dependence on the initialization of the mechanical register. This 
observation is key to the understanding of the working principles of our proposal. 

In order to shed light on this point, we study a single optomechanical system and investigate the dynamical entanglement 
generated when the mechanical (optical) mode is prepared in a thermal (coherent) state. We consider the set of parameters 
listed in Table HJ which are taken from a very recent experiment and represent the current experimental state of the art If2(ill 
and determine the corresponding logarithmic negativity. In Fig. 2 (a) we show the maximum dynamical entanglement achieved 
within the time-window that precedes the reaching of the steady-state conditions corresponding to the choice of parameters 
mentioned above. We clearly see the existence of two distinct regions: in the low temperature region (T<0.02 K), there 
is always an instant of time at which dynamical entanglement is created. In the high temperature region where T>0.02 K, 
dynamical entanglement never arises. In our discussion on discord-enabled entanglement activation, we thus keep such two 
cases distinct. 

Case 1): High temperature activation.- We now consider the case of two independent optomechanical subsystems, whose 
mechanical modes are assumed to be prepared in a thermal state at T > 0.02 K each. In this regime of parameters, the radiation 
pressure interaction is not dynamically entangling and, given the bilocal nature of the overall map, the entanglement in the 
mechanical-vs-optical bipartition has necessarily to be zero. We now show how this scenario changes when initial discord is 
considered. 

Let us take the mechanical modes as prepared in a thermal state at T = 0.4 K each, therefore well in the region of no dynamical 
entanglement. Notice that this is also the temperature of a cryostat. Using the protocol put forward in Ref. [29], which is based 
on the driving of two non-interacting optomechanical cavities with an entangled two-mode squeezed state of light, we can create 
a separable yet non-classically correlated state of the mechanical modes. For instance, using a two-mode drive having squeezing 
parameter r = 1, the joint optomechanical device is such that E = with D = 0.002. We now turn off the driving fields 
and shine the high-intensity uncorrected pumps on each cavity (same parameters as for Fig. 2 (a) but for the temperature) so 
as to check for the possibility to activate entanglement. In Fig. 2 (b) we show how such initial fully mechanical non-classical 
correlations allow to activate optomechanical entanglement in the mirrors-vs-cavity fields bipartition over a time window of the 
order of the decoherence time induced by the Brownian reservoirs. Needless to say, to be faithful to the scheme in Ref. [12], 



TABLE I: Parameters used for the simulations run throughout the manuscript (we assume the same values for each optomechanical subsystem). 
Values taken from Ref. dj]. 



Parameter 


Symbol 


Value 


Mechanical mass 


m 


145 ng 


Mechanical frequency 


w m 


947 KHz 


Cavity length 


L 


25 mm 


Input power 


V 


11 mW 


Cavity-field wavelength 


A 


1064 nm 


Optical damping rate 


k/2tt 


215 KHz 


Mechanical damping rate 


j/2tt 


140 Hz 
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FIG. 3: Low temperature entanglement activation, (a) Optomechanical entanglement created between a single mechanical mode and its 
corresponding cavity field. Each mirror is prepared in a thermal state at the same temperature of its Brownian bath, while the optical mode 
is in a squeezed state. The optomechanical entanglement is shown against time and the squeezing parameter of the optical modes, (b) 
Optomechanical entanglement created in the two mirrors two light modes partition in the following initial conditions: light modes prepared 
in a single-mode squeezed state with squeezing parameter r = 0.5, thermal bath at temperature T — 0.03mK, i) blue line: the mirror are 
prepared in a discorded state; ii) red line: the mirrors are prepared in the factorized state given by the reduced density matrix of the previous 
state. 



we should incorporate the effects of the demon's actions over the register, i.e. the local rotations that are supposed to scramble 
the initial state of the register so as to prevent the activation of entanglement. Formally, this is equivalent to the application of 
the operator ®f_ 1 e' il ' i W^ < +3?i ' to the state of the mechanical oscillators, i.e. the use of the covariance matrix vcjq = RvRF as 

the initial conditions of the dynamical map at the core of our investigation (here R = ©|_i ( ' ' u ' ' 



sin d; cos <&; 



is the symplectic 



transformation associated with a phase-space rotation by the angles di). Experimentally, this is obtained by mismatching the 
time at which the mechanical modes start interacting with the respective pumping field: in the ideal unitary case (where all loss 
mechanisms are neglected) this would result in free dynamics of the mechanical modes for mutually unequal times. We have 
run our simulation by using a locally rotated initial mechanical state and finding that for no value of $i the entanglement in the 
output quadripartite optomechanical state disappears. This is illustrated in the inset of Fig. 2 (b), where we show the output 
entanglement for a sample of 10 3 randomly rotated initial states. We are thus in a position to claim the validity of our Statement 
for Case 1) of our study. 

A remark on the success of the protocol is due. Through the scheme described in Ref. B29I1 and used here to initialize the 
mechanical register, we set quantum correlations in the state of the mechanical oscillators and, at the same time, cool them down 
to effective temperatures lower than that of their Brownian environment. In fact, the mean occupation number of the mechanical 
modes prepared in the separable but discorded state used here is as small as 12, against the 8437 thermal quanta that correspond 
to the initial mechanical state at T = 0.4 K. In order to make sure that the activation of entanglement is guaranteed only by 
the availability of initial discord and not to the effectively cold system we are dealing with, we take a discorded state pr> of 
the mirrors "as warm" as the thermal state at the Brownian bath, i.e. a mechanical state whose single-mode reductions have 
variances of the associated quadratures comparable in magnitude to those of a bipartite thermal state at T = 0.4 K. Formally, 
this leaves reduced single-mode states that are classically squeezed: the associated covariance matrices are fully diagonal with 
only slightly unbalanced non-zero elements which are, nevertheless, larger than the vacuum limit 0.5. From pr> we build up a 
classically correlated bipartite state pec obtained by taking the tensor product of such single-mode reductions (this is equivalent 
to discarding all the correlations in v). We then run the activation protocol using the covariance matrices vo and vcc associated 
with such states finding that only in the former instance the mirrors-vs-optical modes bipartition becomes dynamically entangled 
regardless of the correlations-scrambling local rotations, therefore reinforcing the central role played by discord. 

Case 2): Low temperature activation. - For T < 0.02 K, we know from Fig. 2 (a) that the optomechanical interaction sets 
dynamical entanglement within each optomechanical subsystem. The monogamy arguments used above ensure that also the 
mirrors-vs-cavity modes bipartition is dynamically entangled. This makes the assessment of entanglement activation under 
conditions of low temperature less interesting. Nevertheless, there is still room and significance for our activation procedure. 

In fact, the local dynamical entanglement revealed in Fig. 2 (a) can be fully inhibited when considering, as initial states 
of the optical modes in the activation protocol, two single-mode squeezed states. Fig. 3 (a) shows the optomechanical 
entanglement created as a function of time and the local squeezing parameter for T = 0.03 mK. We notice that for high 
enough values of the local squeezing, the dynamical optomechanical entanglement is completely destroyed. Yet, fully in 
line with the high-temperature case, the presence of discord in the mechanical state leads to dynamical entanglement in the 
mirrors-vs-cavity fields bipartition. To demonstrate this we prepare first the mirrors in a discorded state by using the same 
protocol as before J29tl (again, we use a two-mode driving field with r = 1, although any other choice gives qualitatively 
similar results), run the activation protocol and contrast the results with what is achieved by means of an initial mechanical state 
obtained from the separable discorded one by deleting all the correlations (as before, this delivers a tensor product of classical 
squeezed mechanical states). Fig. 3 (b) shows the output entanglement created in the mirrors-light modes bipartition using 
these initial states. The solid (blue) line is for the entanglement generated using the discorded state, while the dashed (red) one 
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embodies the entanglement generated from the factorized state of the mechanical modes. Clearly, the availability of quantum 
discord activates dynamical entanglement in quite a spectacular way. Such entanglement cannot be inhibited by means of local 
scrambling rotations by the demon, in a way fully in line with the analysis performed previously in the high-temperature case. 
In the long-time limit, finally, the two dynamical conditions converges towards the settlement of steady-state entanglement. The 
validity of the discord-based activation scheme remains thus vindicated. 



Discussion 

We have shown the existence of experimentally accessible conditions under which an entanglement activation protocol requiring 
the pre-availability of quantum correlations in the form of discord holds in the CV scenario as well. While this extends 
significantly the breath of applicability of the framework nicely put forward in Ref. [ 12] for the discrete-variable arena and using 
a rigid unitary gate between the register and the ancillae, we stress that our simulations have been run using the parameters used 
in very recent experiments on optomechanics. This dresses our analysis of further relevance as it embodies a proposal for the 
demonstration of interesting fundamental effects in a mature, realistic and extremely timely experimental setting. 



Methods 

Optomechanical evolution. Here we provide the explicit form of the completely positive dynamical map $, ruling the evolution 
of each optomechanical system. We assume bright input lasers, so that the quantum part of the optomechanical dynamics occurs 
at the level of the fluctuations of the mirrors' and fields' quadratures around their respective classical mean values. The latter 
can be easily determined by means of a standard mean-field approach 13011 . while the former obey quantum Langevin equations 
that can be cast compactly as 



dtU = K&+h l , (i = l,2) 



(2) 



where f^=(8Qi,8Pi,8xi,8yi) is the ordered vector of the fluctuations of the dimensionless quadrature operators Qi = 

(ii\fmiUj % m /h and Pi = pi/^/hmiUJl n for mechanical mode i and Sxi — (Sc\ + 8ci)j\f2, Sjji — i(8c\ — 5ci)/y/2 for the 
corresponding cavity fields. Each 4x4 kernel matrix Kj reads 



K; 
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(i = l,2) 



(3) 



with gi=Xiy/h/ i'^ m i U} m) tne effective coupling rate, c\ + £Aj) the mean amplitude of the i th cavity field and 

the cavity-laser detuning modified by each mirror's mean position. The last term in Eq. (O is the 



A, 



vector of input noise nj(i)=(0, \/2~Ki5x l in (t), ^/2Ki8y\ n {t)), where £i(t) is the zero-mean Langevin force operator ac- 
counting for the Brownian motion affecting the mechanical mode i. For large mechanical quality factors, is correlated as 
(^i(t)^i(t'))=2^lJtBTi5(t—t')/huj l m , with ks the Boltzmann constant and T.- L the temperature of the i th mechanical bath, while 
^in = (^^ii+^^in)/v / 2 and 8y\ n =i{8c} n — 5c\ n )/y/2 are the quadratures of the input noise to a cavity. At typical optical fre- 
quencies, the latter are delta-correlated as {8& in {t)8cll{t l ))=8 ]k S{t - t') with {8&^{t)8c\ n {t'))={8& in {t)8^ n {t'))=Q Vj, k. 
Eqs. (O are solved to get the formal expression for the map $ as 



/ i (i) = $[/ i (0)]=e Kit /i(0) + 



f dt'e Kit 'ni(t-t') 
Jo 



(4) 



Quantum entanglement and discord in two-mode CV states. A general two-mode Gaussian state is fully identified by 

specifying the elements Wy = (xiXj + XjXi) /2 of the covariance matrix v, where we have introduced x = (qi,fti, (72,^2) as 
the vector of the quadrature operators fa and pi (i=l, 2) of the two-mode system. Through Eq. (|4j, the time behavior of the 
four-mode covariance matrix Vf = (fifj + fjfi)/2 of our system can be determined explicitly, thus completely solving the 
dynamical problem. Any covariance matrix can be written as 



T 

7 "2 



(5) 



where oli {0L2) and 7 are 2 x 2 matrices accounting for the local variances of mode 1 (2) and the inter-mode correlations. 
Entanglement can thus be quantified by means of the logarithmic negativity i3"l"l[32tl E= max[0, — In 2(V_)] with i/_ the smallest 
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element of the symplectic spectrum of the partially transposed covariance matrix v p — PvP (with P=1©<7 Z ). The symplectic 
spectrum of a matrix /x is given by the eigenvalues of |i(©? =1 i<7*)/i| ((Tj is the i — x,y,z Pauli matrix). Gaussian discord, on 
the other hand, is calculated as lfl3i[l4ll 

D = f(\fM) - - f(n+) +inf /(V^iti). (6) 

Here, f(x) — (^p-) log[2i=] — (^i^) log[^-i], ^2 = det 02, /x± are the symplectic eigenvalues of v, e = c*i — 7(0:2 + 
t>o) _1 7 T is the Schur complement of a.\, and vq is the covariance matrix of a general single-mode rotated squeezed state, 
over whose parameters inf^., / (y det e) should be evaluated, which can be done analytically for general two-mode covariance 
matrices of Gaussian states IIO. [l4ll . 
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